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Abstract 

Symmetry considerations and a direct, Hubbard-Stratonovich type, deriva- 
tion are used to construct a replica field-theory relevant to the study of the spin 
glass transition of short range models in a magnetic field. A mean-field treat- 
ment reveals that two different types of transitions exist, whenever the replica 
number n is kept larger than zero. The Sherrington-Kirkpatrick critical point 
in zero magnetic field between the paramagnet and replica magnet (a replica 
symmetric phase with a nonzero spin glass order parameter) separates from the 
de Almeida- Thouless line, along which replica symmetry breaking occurs. We 
argue that for studying the de Almeida-Thouless transition around the upper 
critical dimension d = 6, it is necessary to use the generic cubic model with all 
the three bare masses and eight cubic couplings. The critical role n may play is 
also emphasized. To make perturbative calculations feasible, a new representa- 
tion of the cubic interaction is introduced. To illustrate the method, we compute 
the masses in one-loop order. Some technical details and a list of vertex rules 
are presented to help future renormalisation-group calculations. 
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1 Introduction 



The Ising spin glass is the simplest model still incorporating all the complexity 
that more sophisticated disordered systems show up. As such, it has become widely 
studied in the last decades. We focus our attention here to the case where Ising spins 
interact via Gaussian-distributed pair interactions Q . A huge amount of literature 
has accumulated since the seminal paper of Edwards and Anderson nevertheless 
the most important problems — i.e. the nature and complexity of the glassy phase, 
the existence of a transition in nonzero magnetic field, temperature-chaos, etc. — are 
still debated. Consensus has been reached only for mean-field theory, first derived 
considering an infinite number of fully-connected Ising spins || ; its solution by Parisi 
explicitly breaks the replica symmetry, resulting in a picture where the glassy phase 
can be decomposed into a set of infinite number and ultrametrically organised pure 
states 0. Despite all efforts made to go beyond mean-field theory §, which is 
certainly valid in infinite dimensions, finite-dimensional short-ranged systems are 
much less understood. Beside the mean-field picture, an alternative description, the 
so-called "droplet picture" has emerged in a series of papers (a list of them, which 
is certainly not fully complete, is provided as Ref. |5|). This approach claims that 
replica symmetry breaking is an artifact of mean-field theory, and the glassy phase 
consists of only two pure phases related by a global inversion of the spins. The 
droplet theory has gained some support from the field of mathematical physics ||, 
the conclusions, however, remain disputed 0, [||, [|. 

A large amount of numerical workQ has been carried out to resolve the problem, 
finite-size effects and long relaxation times, however, make it difficult to reach a 
definite conclusion. It is clear that analytical methods — especially field-theory, as 
the most powerful of them — are very important to provide reliable results to settle 
this controversy. A direct field-theoretical study of the glassy phase, however, has 
proved to be very hard, due to the complexity of the Gaussian propagators and the 
ubiquity of infrared divergences (see 0] and references therein). A scaling theory 
for the spin glass phase (just below T c and in zero magnetic field) and a proposal 
to handle the infrared problems were put forward in Ref. (!(]], still progress in that 
direction is very slow. 

There is one characteristic of the phase diagram in the mean-field picture which 
is definitely absent in a droplet-like approach, namely the existence of a spin glass 
transition in a uniform magnetic field [jllj , known as the de Almeida-Thouless (AT) 
transition. This question can be studied in the symmetric phase, approaching the 
presumed transition from the high-temperature side, in this way eliminating the 
problems arising from the complexity of the glassy phase. In the language of replica 
field-theory, this will lead to a replica symmetric Lagrangean which is invariant under 
any permutation of the n replicas. 

The prime purpose of this paper is to provide the generic field-theoretical model 
appropriate to a detailed study of the problem raised above, i.e. the existence of 

1 An extensive list of references for numerical simulations in spin glasses can be found in tel. 
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an AT transition. (In a separate publication ||12|| , the crossover region around the 
zero-field critical point is elaborated, and the role that a small magnetic field plays 
in driving the AT transition is investigated.) Here we rediscover, at the mean-field 
level, the importance of the replica number n in the analysis of the AT transition 
fll3| , |j~4|| : for n small but nonvanishing, the AT transition line moves away from 
the zero- field critical point, and an intermediate range of temperature emerges even 
in zero magnetic field. This phase — which can be called, by the extension of the 
concept of Sherrington [15] to continuous n, the replica magnet phase — is replica 
symmetric but still has a lower symmetry than the paramagnetic phase. Hence we 
have two transitions, the first one, in zero field, is an isolated critical point between 
the two replica symmetric phases (paramagnet and replica magnet), whereas the 
second one is a whole line in H — T plane between replica magnet and the replica 
symmetry broken phase. As a result, we can identify the AT transition as the onset 
of instability of the replica magnet phase, and since it has a lower symmetry than 
the paramagnet, we must use a generic replica symmetric Lagrangean to study it by 
field-theoretical methods .0 

Prom Eqs. (|l]) and (Q), one can immediately realize that a perturbational calcula- 
tion based on that Lagrangean is extremely difficult. This is due to the complicated 
interaction term with eight different couplings corresponding to the eight possible cu- 
bic invariants, and also to the non-diagonal Gaussian-propagators with three distinct 
bare masses. To overcome these difficulties, we introduce a new representation of the 
cubic interaction which is associated with a block-diagonalization of the quadratic 
part. The technique proves to be very efficient, as is displayed in our example where 
the one-loop calculation of the mass operator is presented.^ 

The outline of the paper is as follows: In Sec. || the generic cubic Lagrangean for 
the replica symmetric field-theory is set forth, first using only symmetry consider- 
ations. It is then derived, starting from the lattice system, and using the standard 
Hubbard-Stratonovich transformation. Sec. ^ is devoted to an analysis of the zero- 
loop, i.e. mean-field, results. The zero-field transition, first discovered by Sherrington 
and Kirkpatrick Q, proves to be an isolated singularity of the stationary condition, 
with the unique mass vanishing at that point. On the other hand, one of the three 
modes becomes massless along the AT surface, signalling the onset of instability of 
the generic replica symmetric phase. In Sec. |] we define the new set of cubic cou- 
plings. The introduction of this new representation makes it possible to compute 
Feynmann-graphs of a perturbative approach; this is illustrated in a one-loop calcu- 
lation of the three masses in Sec. B. A simple and convenient non-orthogonal basis is 
presented in Appendix A, whereas a detailed list of vertices computed in this basis 
is given in Appendix B. This almost complete table of vertex rules is published here 



2 We keep n > small but finite almost everywhere throughout the paper. This is because we 
want to present formulae for later calculations in the generic cubic model. For this, however, the 
n — > limit is rather tricky, due to the degeneracy of the longitudinal and anomalous modes at zero 
n. 

3 Other methods are also available, like Re plic a Fourier Transform J^J, or the usage of projections 
to the subspaces of the fundamental modes [|L2| , E^ . 
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for later references, an application for an extended renormalisation group study of 
the finite-dimensional AT transition is in progress 18]. 



Cubic replica field-theory for the Ising spin glass in 
nonzero magnetic field 



After the invention of the renormalization group [|T9| , field theoretical representations 
of statistical models, originally defined on a lattice, became a standard way to study 
the behaviour of the systems near phase transitions. The renormalisation group 
made it possible to use a perturbative method, the loop expansion, in low enough 
dimensions, thus providing excellent analytical tools to compute phase diagrams and 
critical properties. The extension to spin glasses came immediately after the replica 
approach had been introduced by Edwards and Anderson ||, transforming the orig- 
inally inhomogeneous system into a homogeneous one. A Ginzburg-Landau- Wilson 



continuum model was first proposed 25 years ago [20, 21], then further investigated 
by renormalisation methods p2| , |23| . Its cubic Lagrangean was derived — for the 
Ising case — by Bray and Moore [£4[ via the Hubbard-Stratonovich transformation. 
Two of the present authors applied the same field theoretical model in their effort 
to go beyond the mean-field results, and understand the glassy phase of the finite- 
dimensional short range Ising spin glass The magnetic field was always zero 



in the above works, with the only exception of [21], where it was introduced by a cou- 
pling to the magnetization, leaving the Lagrangean unchanged for the part relevant 
to the spin glass transition. 

Field-theoretical models can be constructed by means of symmetry arguments, 
building up the Lagrangean from all the possible invariants of the relevant symmetry 
group of the system. For an Ising spin glass, the fields depend on a pair of replicas, 
(j) ! 3 = (f>P a with 4> aa = 0, and, as a consequence of the replica trick, any permutation 
of the n replicas leaves the Lagrangean unchanged. Discarding all the terms higher in 
the order of the 0's than cubic, we arrive at the following generic replica symmetric 
Lagrangean after a search of all the quadratic and cubic invariants: 

£ = £ (2) +£ (3) 



where 

c® = ly 

2 ^ 
p 



i( W ) 2 + m 1 )^f^ p + m 2 £, 



a(3 



aft"/ 



>p 7 ^ 7 p 



+ m 3 £ 

a/375 



(1) 
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and 



P1P2P3 



£ (3) = -^ £' -iE«c+- 2 E<«« 

v , , _ , , _ , , . _ a/3-y a/3 



+«*£ «C + ^ £ + 



E 

a/37 aftjS a/376 



+^E<c<+- 7 E E «c 

a/375 a/375^ a/375/iv 



(2) 



Momentum summations in the above formulae are over the reciprocal vectors of a 
d-dimensional hypercubic lattice with lattice spacing a, consisting of infinitely many 
sites N in the thermodynamic limit. The prime in Eq. (||) means the constraint of 
momentum conservation, pi + P2 + P3 = 0. Neglecting the fluctuations of fields with 
wavelength much smaller than the range pa of the exchange interaction between 
the spins, we confine the relevant values of momentums in Eqs. (|l],||) around the 
center of the Brillouin zone: p < A/ pa. The momentum cutoff A <C 1 makes it 
possible to expand the nonlocal quadratic coupling in C^ 2 \ and, as it is common in 
field-theoretical studies of phase transitions, we stop after the first two terms. (The 
coupling constant of the kinetic term in Eq. (|l|) can be set equal to \ without loss 
of generality. See also later.) 

In zero magnetic field, and in the high-temperature phase where the spin glass 
order parameter is zero, all the couplings but mi and w\ disappear. In this section, 
we want to find out the order parameter dependence of the couplings defining our 
Lagrangean. We are especially interested in the general form of the replica field- 
theory suitable for studying the de Almeida-Thouless transition in finite dimensions. 
Our starting point is a standard Edwards- Anderson-like g] model for N Ising spins 
on a d-dimensional hypercubic lattice, with a long but finite-ranged interaction: 

w = (3) 

(y) V 

The notation /y = / ^ r ' pa rj ^ was introduced in the above equation, with the smooth 

positive function f(x) which takes the value 1 for x < 1, and decays to zero suffi- 
ciently fast for x > 1, thereby cutting off the interaction around |rj — | ~ pa. 
Here z = p d is effectively the coordination number, i.e. the number of spins within 
the interaction radius; expanding quantities in terms of its negative powers will gen- 
erate the loop-expansion in the replica field-theory. Jij are independent, Gaussian 
distributed random variables with mean zero and variance A 2 , and a homogeneous 
magnetic field H was also included. Summations are over the pairs (ij) of lattice 
sites in the first sum, while over the lattice sites in the second one. 



4 A replica symmetric treatment of the ordered phase was carried out in Refs. JiJ], ^3. In this 
case, the finite order parameter gives rise to a quadratic Lagrangean, even in zero field, which is a 



special case of £' 2 '. 
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In the spirit of the replica trick, we want to compute quantities like the averaged 
replicated partition function Z n for some positive integer n, finally deducing spin 
glass behaviour from the n — > continuum limit. Averaging first over the quenched 
disorder results in an effective replica Hamiltonian depending on the spins Sf , a = 



(a/3) y 



1 / A A 2 j?2 



y = - ^^p) /?.. A Hubbard-Stratonovich transformation can help us to get rid of 
the four-spin interaction term; the price we have to pay for that is the introduction 
of integrals over the "fields" Qf 13 '- 





" n / 


exp ^ 









(4) 



(a/3) 



The boldfaced vector and matrix notations in Eq. (Q) for and K _1 , respectively, 
are to simplify the nonlocal term in the formula, whereas the one-spin partition 
function is defined as follows: 



H 



{5 1 V/3) 



(5) 



To construct a field-theory appropriate for a perturbation expansion around the 
mean-field solution, i.e. around the infinite range model p — > oo, we separate 
into its homogeneous, non-fluctuating (mean-field) part, and into its fluctuating part: 

Qf = Q°>P + <f>f. ( 6 ) 

When expressed in terms of the $>'s, the exponent in Eq.(|^) will be called — C, 
and it can be expanded up to any desired order. Turning to a more convenient 
representation of the fields in momentum space, contributions up to cubic order 
have the following forms: 



£(0) = N 



V^(Q^) 2 -ln C 



(a/3) 



(a/3) 



e -lQa/3 _ (( S a SP)) 



°p=0 > 



£(2) = ^E E <t>fMaw(pap)<t>Z S p 



(7) 

(8) 
(9) 



P (a/3),( 7 5) 



5 Throughout the paper, we use different notations for summations over distinct pairs, — 

(a/3) a</3 

and unrestricted sums, XI ~ S 5Z- 

a/3 a (3 
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and finally 



£(3) = "£7f ^ S ^C" ( 10 ) 

V P1P2P3 (a/3),(yS),(fj,v) 

A Boltzmann-weight with Q a ^, instead of , is understood in the definitions of 
( and the one-site effective average ((■■■)) in Eqs. (@^|). _1 is essentially the 
temperature squared, or more precisely: 

a \ 2 (• 



e= {kf) j^ d r - ( U ) 

The momentum- dependent mass, and the momentum- independent cubic coupling 
operators are defined as follows: 

M a ^ s (pap) = K-Xp, r s ~ [((S a S?S^S s )) - ((S<*SP))((£nS 5 ))] , (12) 
Waprfw = {{S a S p S' r S s S f "S u )) - {{S a sP)){{S"<S s S»S u )) 

- ((S^S s ))((S a S p S^S u )) - {{S' 1 S u )){{S a S /3 S"'S 5 )) (13) 
+ 2((S a S l3 ))((S^S s }}((S^ l S u )}. 

The Kronecker delta in Eq. (^) represents the n(n — l)/2-dimensional unit matrix, 
whose prefactor comes from the Fourier-transform of Kif. 



ij 



The arrow in the above formula means the double limiting procedure of the thermo- 
dynamic limit (N — > oo), and continuum limit (a — > 0). The theory resulting then is 
a field-theory with all the thermodynamic functions scaling correctly with N, and the 
lattice constant a disappearing from the momentum integrals after introducing the 
new variable p = pap, with the upper momentum-cutoff becoming A in p space. The 
range p of the original interaction, however, survives: a perturbative loop-expansion 
can be generated where every loop in a Feynmann-diagram contributes a z^ 1 = p 
factor. 

It is rather common in field-theoretical studies to normalize the fields such that 
the kinetic term in the Gaussian part of the Lagrangean (that proportional to the 
squared momentum) be exactly p 2 times the unit matrix. This can be simply reached 
after expanding K~ l , Eq. (14), and introducing the new fields by 



y p< 



(15) 
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where 



A corresponding redefinition of the mass operator and cubic interaction, 

\m^M and \]V^W, (17) 

cr c? 

leaves the form of Eqs. (p| |To|) unchanged. Neglecting short-wavelength fluctuations, 
the mass-operator in Eq. (|12|) can be expanded for p <C 1. The commonly used 
truncation at the kinetic term provides: 

M a ^ 5 (p) = C{d) [sf MS - e{((S a SCS^S 5 )) - ((s^s^as-rs 6 ))) 

+ P 2 fe, (18) 

where C(d) = 2d J f(r) 2 d d r/ J r 2 f(r) 2 d d r is a smooth function of dimensionality, 
but independent of the temperature and magnetic field. As such, its concrete value 
is irrelevant, and a simple adjustment of the cutoff function f(r) can make it equal 
to unity. 

A replica symmetric field-theory — for the study of the massive high-temperature 
phase, and/or the massless critical manifolds — can be obtained by choosing a replica 
symmetric mean-field value Q a p = Q in Eq. (pi). The stationarity condition C^ 1 ' = 
gives us an implicit equation for Q (see Eq. (ra)): 

Tr ^"^e^ic/?) O^^+h 7 S a s° j 
e- 1 Q = ((S a SP)) = — r — —g r , (19) 

{s a } V / 

Q enters the definition of the mass operator, Eq. (^), and the cubic interaction 
operator, Eq. (|l3|), through the Boltzmann- weight in the averages ((...)). Replica 
symmetry is induced also for these operators, resulting in the three different compo- 
nents of the mass: 

M a p, a p(p) = M 1 +p 2 , 

M Q7i/37 (p) = M 2 , (20) 
M Q/3 , 75 (p) = M 3 ; 

and the eight different components of the cubic interaction operator: 

Wapfi^jat = Wi , W a p :a p :a p = W 2 , W a p )a p )0n = W 3 , W a p, a p }1 S = W4 , 
W a p ta% ps = W 5 , W a p,orf,aS = W 6 , Way^Sfi = Wj , W a p, 7 S,ixu = W 8 ■ (21) 
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= Wi- 


3W 5 + 3W 7 - W 8 




U>2 


1 

= -w 2 


- 3W 3 + -W 4 + 3W 5 + 2W 6 


- 6W 7 + 2W 8 


w 3 


= 3W 3 - 


- 3W 4 - 6W5 - 3W 6 + 15W 7 


-6W 8 


W4 


= \w. 


- \w 7 + ^W 8 




W5 


= 3W 5 - 


- 6W 7 + 3W 8 




w (i 


= w 6 - 


3W 7 + 2W 8 




W 7 


= -w 7 

2 


--W 8 
2 






1 

= 0^8 







Table 1: The relationship between the cubic couplings w's of Eq. (|2|) and W's of 
Eq. m. 



of Eq. @, together with Eq. (|20|), is equivalent with that of Eq. (jlj), provided 
the two sets of masses are related by the following expressions: 

mi = ^(M 1 -2M 2 + M 3 ), (22) 
m 2 =M 2 -M 3 , (23) 
m 3 = -M 3 . (24) 

Similarly, Eqs. (||) and ( |l~0| ) are two different representations of C^ 3 \ Using Eq. (21), 
a one to one correspondence between the two sets of cubic couplings, w's and W's, 
can be deduced by a somewhat lengthy but elementary calculation. The results are 
listed in Table 1. 

3 Analysis of the stationarity conditions and bare masses 



It is easy to recognize that, after a simple rescaling of the temperature, Eq. (|19| ) coin- 
cides with the replica symmetric mean-field equation of Sherrington and Kirkpatrick 
(SK) |3|, ^] for the order parameter of the Ising spin glass on a fully-connected lat- 
tice. This may not be surprising: the most direct way to define mean-field theory on 
a d-dimensional hypercubic lattice is letting p, the range of interaction, go to infinity, 
thus neglecting all the loop corrections to the equation of state |^6|, ^?J. The solu- 
tion Q of Eq. ( |l^ ) has, however, an application that goes beyond mean-field theory: 
it enters the mass operator and interaction components in the formulae Eqs. (|i~8| ) 
and (|l~3l), respectively, through the effective average ((...)). Although a field-theory 
emerging from this procedure has a direct connection to the original parameters, 
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such as temperature, magnetic field and also replica number n, renormalisation will 
reshuffle the masses and couplings, possibly forcing them into some fixed point. 
Following Refs. P, 



, Eq. (19) can be cast into a more convenient form: 



J 4=e-^r tanh 2 (^/Qu + H/kT) cosh n (^/Qu + H/kT) 



J ^e-ir cosh n (^Qn + H/kT) 



(25) 



tank 2 (^/Qu + H/kT). 



The shorthand notation 



j jhi e -^ tanh *(. _ _ ) cosh n (V^w + H/kT) 



tanh%..) 



j du e -- cosh ™(yQu + H/kT) 



(26) 



has been introduced for later use. By Eq. (|25|), Q is implicitly given as a function 
of temperature, magnetic field and replica number. One can easily find the SK spin 
glass transition point (0 = 1, H = 0) as an isolated singularity for any given n close 
to zero. (Keeping n finite is for later use. At the moment, we must notice that this 
singularity is rather unaffected by the n — > limit.) The relevant, positive, solution 
for H = is 



Q 



o 



2-n 



t + 



n—3 



(n-2) 2 + 3(n-2) a 



t Z +... 



for t > 0, 
for t < 0; 



(27) 



where the new temperature scale t = 0" 1 — 1 > (< 0) in the disordered (spin 
glass) phase, respectively. As displayed in Figure 1, one can join up smoothly the 
two regimes by-passing the critical point, like in ordinary critical phenomena.^ 

To go beyond a mean-field solution, and build loops, one must have well-defined 
Gaussian propagators, i.e. the eigenvalues of the mass operator must be non-negative. 



A generic replica symmetric mass operator, like that in Eq. (20), was diagonalised 
years ago |23|, 24 1 with the following expressions for the eigenvalues (which will play 



the role of bare masses here): 



r R = M x - 2M 2 + M 3 , 
r A = Mi + (n - 4)M 2 - 

r L = Mi + 2(n - 2)M 2 + 



[n 



3)M 3 



(n-2)(n-3) 



M,. 



(28) 
(29) 

(30) 



The indices R, A and L stand for replicon, anomalous and longitudinal, respectively, 



each referring to the corresponding family of eigenmodes (see Ref . [ 24 ] and also later 

6 It is obvious from Eq. (|25|) that Q = is always a solution for zero magnetic field, independently 
of t and n. The Q — solution for t < 0, however, defines the negative Q branch (starting from 
point G in Figure 1(a) with Q — 0, and following the dotted curve, you end up at P with Q < 0) 
which is non-physical. The two branches meet at the SK transition point. 
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(a) 



G P 

(b) 



Figure 1: (a): Mean- field phase diagram for n = 0. Starting from the paramagnetic 
state P with Q = 0, and following the dotted curve, the SK critical point (t = H = 0) 
can be by-passed ending in the glassy state G with Q > 0. (b): Mean- field phase 
diagram for n > 0. There is a temperature range even in zero field, *at < t < 0, 
where Q is positive and the replica symmetric state is stable. 



sections). From Eqs. ( jig ) and (20), and using the stationarity condition (|19|), the 
following expressions are obtained for the bare masses: 



t 



+ 2Q 



1 



l+t " 1+t 
r A = r R + (n - 2 
n — 1 



tanh^i^/Qu + H/kT), 



ri = r R + 



—— tanh 4 ( v / Qu + H/kT) - Q , 

± ~\- t 

(n - 4)—— tanh 4 ( v / gu + H/kT) - 4Q + n(l + t)Q 2 
1 H~ c 



An expansion below the transition provides, in zero field, 

r R = 2(1+ 2 +...)(-*) + (-3 + --0* +•• 
r A = (l + ^ + ...)(-t) + (--^ + ...)t 2 + ... 



rL = (-0 + (-^ + ---)i 2 + 



(31) 
(32) 

• (33) 

(34) 
(35) 
(36) 



(The n-dependent coefficients are displayed in an expanded form for showing clearly 
the signs for small n. The complete n-dependence, however, is easily found.) The 
hitherto degenerate masses split when passing the SK transition singularity (where 
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they are zero) and emerge positively in the spin glass phase for any n > 0. 
starts, however, with a small slope proportional to n, and becomes massless again 
at the AT surface t A T = — §ra + • • • , where instability of the replica symmetric phase 
begins. This result was first derived by Kondor jL3|, the aspect we wish to emphasize 
here is the existence of an intermediary temperature range where replica symmetry 
persists, see Figure 1, though, as a result of a nonzero Q, the level of symmetry is 
lower, leading, at the mean-field level, to the splitting of the bare masses. 

By definition, m = on the AT-surface. From Eqs. ( p5[ ) and (|3l"|), the magnetic 
field can be expressed as a double series for small —t and n; the leading, cubic, term 
is as follows: 

(H/kT) 2 = h-tf - -nt 2 = ~ln 2 (t - t AT ) + \n{t - t AT ) 2 - \{t - t AT f. 
o 4 8 2 o 

This can be cast into a scaling form 

(H/kT) 2 = t T ^), i,n-0, (37) 

where the exponents have now their mean-field values r = 3 and k = 1, and cp(. . . ) 
is the scaling function characterising the AT-surface. 

On the basis of the above mean-field analysis, we can define two different, though 
both replica symmetric, cubic field-theories relevant to describe spin glass transitions 
of different types in low enough dimensions (d certainly smaller than eight): 

• The zero-field model — H = 0, t > and n > — implies, through Eq. (|27j), 
Q = 0, leading to degenerate bare masses tr = ta = ri, see Eqs. (31), ( |32| ) 



and (pi). It is easy to verify using Eqs. (Ill) and pi) that all the W's but W x 



are zero, and the same is true, by Table 1, for the small w's. Let us put down 
explicitly the definition of the cubic field theory for the zero-field transition 
from the paramagnetic phase: 

r =rr? = r a = rr; 

(38) 

w =wi, Wi = 0, i = 2, . . . ,8. 

At the mean-field level, a massless state is reached along the critical line H = 0, 
t = and n > 0; the replica number being a rather innocent parameter around 
zero. We expect this property to remain for short-ranged (finite-dimensional) 
systems, and indeed, e-expansion results have supported this idea Jp8| . (This 
kind of replica field-theory was studied in all the Refs. gl], ||, || 




The second model, which is in fact the most general cubic field-theory with an 
unbroken replica symmetry, has all the three masses mi and eight couplings Wi 
different. At the mean-field level, it corresponds to a nonzero Q which is always 
such when a magnetic field is switched on. More surprisingly, however, there 
is a whole range of temperatures £at < t < even in zero field where the bare 
parameters correspond to this more general model. Criticality is induced, at 
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least at the mean-field level, by the masslessness of the replicon mode on the de 
Almeida-Thouless surface. The replica number n is now a crucially important 
parameter; a fact clearly shown by the scaling formula (|37]). How fluctuations 
will modify this picture is a prime problem in spin glass theory. There has been 
two attempts, at least to our knowledge, adressing this question (l?], p9{| . In 
Ref. [17] fluctuations were restricted to the replicon subspace; in our language 
this means for the masses that = r*£ = oo and tr critical, while all the cubic 
couplings were zero but w\ and W2- The effect of a small magnetic field was 
introduced in Ref. |2£| by shifting the bare masses from their zero-field values 
(more precisely, beside mi, rri2 became massive too), the couplings remained, 
however, unchanged. The replica number was effectively set to zero in these 
works. Here we wish to emphasize the role n may play in a search for a de 
Almeida-Thouless transition in finite-dimensional systems. 



4 A canonical representation of the cubic interaction 

In a field-theory with more than one mass, as in the generic replica symmetric 
system introduced in the previous sections, critical manifolds can be classified by 
their massless eigenmodes. In this case, it is more convenient to use directly the 
eigenvalues of the mass operator (t\r, ta and r^; Eqs. (28), (2£) and (30)), instead 
of the sets m's or M's. At that point it is natural to ask what will happen with the 
interaction vertices whose legs join the, by this time block-diagonalized, propagators. 
We show in this section how the transformation that block-diagonalizes the quadratic 
part of the Lagrangean into "modes" induces a new set of cubic couplings describing 
how these modes interact. Our replicated field-theory becomes more tractable after 
using these "canonical" cubic parameters. 

The \n{n — l)-dimensional vectorspace spanned by the two-replica fields (f) a/3 
has the simple structure being a direct sum of the subspaces called longitudinal, 
anomalous and replicon. Their definitions are as follows (<j) al3 = (j^ a and cf) aa = 
are understood everywhere, of course): 



The longitudinal (L) subspace consists of replica symmetric vectors, i.e. inde- 
pendent of replica indices. Each element from this subspace corresponds to a 
scalar (f) : 

4>f = <t>, (39) 
and it is obviously one-dimensional. 

Any element of the anomalous (A) subspace can be represented by a one-replica 
field 4> a , i.e. by a vector restricted, however, by the condition 



5> 



0. (40) 
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A generic anomalous field can now be written as 

<^t = \{^ + ^)- (41) 

As a result of condition (|4^), the anomalous subspace is n — 1-dimensional. 

• True two-replica fields, loosely speaking tensors, constitute the replicon (R) 
subspace with the restriction 

¥r = for any a = 1 . . . n. (42) 

Prom the n equations above follows that the number of independent (ff^ 
is \n{n — 1) — n = \n(n — 3), rendering the replicon subspace ^n(n — 3)- 
dimensional. 

A generic field 4> a/3 can always be decomposed into the sum 

= if + <t>°£ + (43) 

It is straightforward to see that the subspaces defined above give the exact di- 
agonalisation of a generic replica symmetric matrix, as defined in the equations of 
©, i-e. 

M ^a^t = ntf , i = L,A,R , (44) 

(r/S) 

the eigenvalues given in Eqs. (|30[), ( |29| ) and (28), respectively. For a generalization 
to higher order operators, the matrix-element of M between two arbitrary vectors 
(j) al3 and tp al3 will be computed after having represented them by the longitudinal 
scalar, anomalous vector and replicon tensor, as explained above (see Eqs. (|3^), (|4l|) 
and fl45|)). An expression in terms of the three second-order invariants <fnp, Y^a 4' a ' t P a 
and Eq/3 4>riP'r arises: 

(a/3). (7*) , . 

(45) 



1 a/3, a/3 , n — 2 ST^ ,a , a , n(n — 1 ) 

2 r «L^A +— r ^L^ ^ + — o — rL ^ 



a/3 



(to get the anomalous term, the restriction (40) has been used). Except for numerical 
factors, the corresponding eigenvalues appear as the coefficients of the three possible 
second order invariants; viz. RR, AA and LL. 

The most important point we can learn from Eq. (|45| ) is a complete factorization 
of a matrixelement, provided 4> and ip are chosen from one of the subspaces L, A or 
R; i.e. 



M {<t>M 



{expression of mass components} x {invariant composed of <j) and tp}. (46) 
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Generalization to a cubic replica symmetric operator, as defined in Eq. is 
straightforward. The analogue of the matrixelement can be easily defined by 

w m,x}= E W^s^r 6 ^- (47) 

Taking each of the fields 4>, ip and x from one of the subspaces L, A or R, the nonzero 
values obtained can be listed as follows: 



Wrrr =9iJ2 a/ V V° + \g2 E ^ 

a/37 a p 

W RRA = g 3 5> a/ W, 



a/3 a/3 
A 1 



a/3 



W RRL =g±^<f> a0 ^X, 

a/3 

W RAA =g 5 J2<t> a W, (48) 

a/3 

W AAA = g 6 ^2^ a X°, 

a 

W AAL = g 7 Y,^ a X, 

a 

W LLL = g$ (/npx- 

Symmetry makes the remaining W RA l, W R ll and W A ll all identically zero. As for 
the masses, a complete factorization occurs in the above formula, except the RRR 
vertex.[] We prefer the name "canonical" for the set of cubic parameters g emerging 
in the above formulas, as they are, in some sense, an extension of the notion of 
eigenvalues to the cubic interaction term. After a somewhat lengthy calculation, 
we obtained the set of equations for the gi in terms of the Wj, i = 1, . . . ,8, which 
are the counterparts of Eqs. fl2S|), (|29|) and (|30|). (We omit to display these rather 
complicated, and not very instructive, expressions here; they can be easily obtained 
from Eqs. (49a-h) below and using Table 1.) 

Comparing Eqs. (|2|) and a one to one correspondence between a wi and a 
gi, i = 1, . . . , 8, is obvious. When expressing the g's in terms of the w's, instead of 
the W'b, not only the formulas become simpler, but a clear R — ► A — » L hierarchy 
emerges: 

gi =w 1 , (49a) 

92 = 2w 2 , (49b) 

n — 2 

93 = -wi +w 2 -\ — w 3 , (49c) 



7 The reason for that is the two different cubic invariants we can construct from replicon fields; 
see the first line of Eq. (fig}) . 
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n — 1 n(n — 1) 

g A = - Wl + W2 -\ — -| w 4 , 

n-4 1 n-2 (re - 2) 2 
95 = — : — wi + -w 2 H — w 3 H — w 5 , 



6 



3n-8 



.96 



97 



4 

n-2 



n-4 



-u>2 + 
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(n-2)(re-4) (n-2) 5 



(49d) 
(49e) 



1^5 + 



(n - 2) s 



n — 4 2n — 3 n(n — 1) 
— - — wi + w 2 -\ — u>3 H ^ W4 



(n-2)(2n-3) n-l)(n-2) n(n-l)(n-2) 
6 2 6 

gs = n(n — 1) (n — 2)^! + w 2 + (n — 1)^3 + n(n — l)u) 4 
+ (n — l) 2 u>5 + (n — 1) 2 ^6 + n(n — l) 2 u>7 + n 2 (n — l) 2 tt>8 



(49f) 



(49^ 



(49h) 



5 Illustration of the technique: one-loop calculation of 
the masses 

We apply the standard definition of the mass operator, i.e. it is the zero momentum 
limit of the inverse of the two-point function: 



r a /3 i7< 5 = hm [G ^P)]^ 



(50) 



(51) 



where the connected two-point function G is the average 

taken with the statistical weig ; see Eqs. (|]) and @. Dyson's 

equation for T allows us to compute it perturbatively: 

r a /3, 7 5 = M a/3aS (p = 0) - S Q/ g i75 (p = 0), (52) 

where the mean- field mass operator M a p^§ has been defined in Eqs. (^) and (pCf), 
whereas the self-energy £ contains all the one-particle irreducible graphs to the two- 
point function, with external lines omitted. Up to leading, one-loop, order it is given 
as the simple "bubble" diagram: 



S a /3, 7 5(p = 0) 



1 



d d p 



^ W a p,a'P',a''l3''G°,^ tl , s ,(p)G^ ll3 „ tl ,, s „(p)WjS,j'5',j''S'' (53) 



2z J (27r) d 
where the free propagator G^ ' is defined by 



(0) 



G 



(0) 



(54) 
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see Eqs. ( fig ) and (|20|), and the bare vertices, W's, have been introduced in Eqs. 
(H), (|D and (U). 

To compute the replica sum in Eq. ( 53 ) , we must overcome the problem of having 
non-diagonal free propagators. This project can be easily accomplished by calculat- 
ing Yj m n: instead of Yj a ^^g, defined as: 

S w>n = 0m S a/J, 7 «J 4>n, (55) 

(q/3),( 7 5) 

4> m and 4> n taken from the non-orthogonal basis discussed in Appendix A. Exploit- 
ing the completeness of both this basis and its biorthogonal counterpart^, we can 
transform S mj „ into a representation with diagonal free propagators: 

Gffii^ £ rjG^ s ^ = G^5% n (56) 

(a/3), (7-5) 

(a "tilde" refers always to a member of the reciprocal basis). Gm > as it is the 
eigenvalue of the free propagator matrix, can be simply related, through Eq. (54), 
to one of the three bare masses of Eqs. (|28|), (|29|) and (j3C|): 

G$(p) = — (57) 

fm = rR> or vl depending on the subspace m belongs to (see Appendix A). 

We can now propose a simple graphical representation for S mj „ by introducing 
arrowed lines for the free propagators joining interaction vertices W: 

n(°) <<=> m 

w 

Using the convention of Appendix B concerning the meaning of inward and outward 
arrows, we can draw for E m;n : 




which can be spelled out explicitly as 

r-A A d~ 

2z 



rj/ A ?Sd E W m , ml)W/ W n>m ^^ G^{p)G^,(p). (58) 



3 What we use in the derivation of this formula is the decomposition of the unit operator: 5 a a lT a 



Em C ft = J2r, 
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(Summations are over the set of n(n — l)/2 modes of Appendix A.) Orthogonality of 
the different subspaces restricts the number of nonzero elements of the matrix S m n 
to the cases where m and n belong to the same family R, A or L. To compute the 
three eigenvalues of the self-energy, we can make the simplest possible choices for m 
and n, i.e. 

v v aP-P (n- l)(n-2) 2 (n-3) 

=^R #G«0 #(/«/) = 4 hR ' ^9) 

(a/3) 

v v \^a<*Pa<*P n(n-l)(n-2) 

E(aO,(aO = S A ^ ^)%) = 4 , (60) 

(a/3) 

%),(L) =^lJ2 tfvtfL) = n(n 2 ^ Sl • ( 61 ) 

(a/3) 

(The computation of the scalar products above is relatively easy using the defini- 
tions of the basis functions in Appendix A.) An extensive use of the table of cubic 
vertices W mm i m" m Appendix B makes it possible to compute the left-hand sides; 
the feasibility of the calculation is, however, due to the selection rule we explain in 
that appendix. The results can be summarized by displaying the eigenmodes of the 
mass operator T, by means of Eq. (52), valid to first order in 1/z: 



Tr = r R 



re 4 - 8n 3 + 19n 2 - An - 16 , 2(3n 2 - 15n + 16) 

9l + —, 777 t^t- 9i92 



(n-l)(n-2) 2 JL (n-l)(n-2) 2 



n 6 - 9n z + 26n - 22 2 
+ 2(n-l)(n-2) 2 92 



Irr 



8(n-l)(n- 4) 2 8 2 16 2 \ 

+ n(n-2) 2 53 lRA + T^T) 54 ^ + J^W 95 J ' (62) 
f2(n-3)(n-4) 2r 16n(n - 3) 2r 
FA = rA " \ („ - 2) 2 53 /flR + (n-l)(n-2) 2 95 IrA 

+ ^Bif 91 ^ + n(n - l)(n - 2) 2 9 " IaL ) 5 (63) 
( 2(n - 3) 2 16 2r 4 2 \ 

Fi = ri " 94 + <^ * lAA + n 3 (n-l) 3 58 7 " j • (64) 

To help the reader to understand the structure of the corrections to the masses, we 
introduced the short-cut notation 



1 f A d d p 1 1 



z J (2-7r) d r s + p 2 ry + p 2 



for the momentum integrals; s and s' correspond to one of the subspaces R, A or L. 
By means of the transformation rules between the g and w couplings, Eqs. (49a-h), 
the above equations for the masses can be easily expressed in terms of the u>'s too. 
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Appendix 

A A simple non-orthogonal basis 

For applying the canonical vertices defined in Eqs. (47) and (|48|), it is necessary 
to introduce a basis in each of the subspaces; an obviously nonunique task. The 
non-orthogonal basis defined below is not only the simplest^], but cubic vertices 
evaluated in this basis will have a remarkable property, a kind of a selection rule 
involving replica numbers, making computation of Feynmann-graphs feasible (see 
Appendix B). 

A member of this non-orthogonal basis will be denoted by (thn , whereas its 
biorthogonalf ! counterpart, a member of the reciprocal basis, as , where m stands 
for the modes in the subspaces L, A and R as follows: 

• The L subspace is one-dimensional, i.e. m <-> (L). 

• m runs the single replica numbers, except one (which we choose the n th ), for 
the n — 1-dimensional A space: m «-> (/i), \i = 1, . . . ,n — 1. 

• In case of replicon modes, m corresponds to a pair of replicas, m <-» (pis), with 
[i,v = 1, . . . , n — 1 and \i ^ v. To ensure the correct dimensionality n(n — 3)/2 
imposed by condition (f42j), we have to pick out an arbitrarily chosen pair (flu), 
giving for the number of replicon modes: 

(n — l)(n — 2) n(n — 3) 

2 ~ 2 ' 

To sum up, there are two types of replicon modes: 

\(^u), n,v = 1,... ,n- 1; H,v^fl or v, 

m <-> < _ _ (65) 

l(/i/i) or (^v), /i = 1, . . . , n — 1; P^i^l 1 or v. 

In what follows we want to collect the results, omitting any proof. 
L subspace: 

2 

0(5 = hL) = i; 4>°l) = 4>(l) = n{jl _ Vj - ( 66 ) 



9 This time, simplicity and orthogonality contradict each other. An orthogonal, still rather com- 
plicated system was proposed years ago BQ |. 

10 Biorthogonality has the usual definition J2(a/3) 't'^'K^i = ^mm'- 
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A subspace: 



1 if a, ft 7^ fj,, 



n-1 



■2f2 if a = M or /? = M ; ^ ™ 2 (™ " 2 ) TO ' ^ 
The one-replica objects, introduced in Eq. (pl|), representing them are 

if a 7^ /i, n, 



I if a ^ /x, 

(n — 1) if a = /i; 



/ 4 \ if a = u, 

n(n— 2) 

! if a = n. 



„ n(n— 2) 



R subspace: 



(67) 



(68) 



V) 



ra—3 
2 



if Q!, /3 7^ fl, V , 

if there is one common replica index of the 
two pairs a, [3 and fi, v , 

k ( n 2 K n 3 ) ^ ^ wo p a i r s Q!) ^ a nd ^ j, are identical. 



(69) 



In the reciprocal basis, we have different forms for the two types of replicon modes 
defined in Eq. (65): 



type-I 



type-II 



V) 



(n-lj(n-2) if = (^) OT OT 

- (ra _i) 2 (ra _ 2) if («/3) = (m?) or (/H or ( vn )i 
otherwise; 



(mm) 



(n-l)(n-2) 
2 



if (a/3) = (/u/x) or (zm), 



~ (n-l)(n-2) if ( a #) = W OT (^ n )' 

otherwise. 



(70) 



(71) 



B Vertex rules 



For a cubic vertex with 



ip = m / and x = 0m" in Eq. 



(to, m' and to" 



referring to the modes introduced in Appendix A), a simple graphical representation 
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can be given, namely 



inward (outward) arrows correspond to ordinary (reciprocal) basis functions, respec- 
tively. Such vertices have the remarkable property, a kind of a selection rule, that 
the replica numbers attached to the mode m"P| must occur either in m or in m!\ 
otherwise the vertex is zero. 

Hereinafter we give a list of the nonzero vertices. To confine the extent of the 
paper, vertices with replicon modes of type-II will also be omitted, although they are 
available; these vertices are necessary only for a calculation higher order than one- 
loop. The presentation follows the order introduced in Eq. (f48|), different symbols 
are used to indicate different replicas. 



RRR,: 



x n—1. . 



n 



■[(n-l)(n-4)si+8j] 



1 



2(n - 2) 



2(n 



_L_ [( n 2 _ 9n + U ) gi _ ( n 2 _ 6n + 7 ) 52 ] 



[2(3n 2 - 15n + 16)51 + (n 3 - 9n 2 + 26n - 22)g 2 ] 



11 Tt is a single number (^1) if m" is an anomalous mode, whereas replicon modes are labeled by a 
pair of replicas, as explained in Appendix A, (fiv) or (fift). There is, of course, no restriction if m" 
is the longitudinal mode. 
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RRA: 



RAA: 



m 2(n - If 



n(n - 2) 



53 



(n- l) 2 (n- 4) 
n(n - 2) 



■53 



t _ (n- l) 2 (n- 3) (n- 4) 
n(n - 2) 



93 



= %3 



n 



n-2 



.93 



(n-l)(n-4) 
n-2 



53 



RRL: 

± 2(n - 2) 



n 



54 



l (n-2)(n-3) 



n 



■54 



l (n-2) 2 (n-3) 



n 



54 



254 



n-2 



55 



m_ 2(n- 1) 
n-2 



■55 



NTjt 2(n-l)(n-3) 
^ _ ^= Z n 55 



n-2 



m_ 2(n - 3) 
n-2 



55 
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2n 2 

95 



(n-l)(n-2) 

AAA: 

n- 2 96 /^=~^ 

AAL: 

9. 

v 97 AT^= 2 97 

n — 1 



4 



LLL: 



n(n — 1) 
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